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Molecular Rotations in Matter-Wave Interferometry 

Benjamin A. Sticklers’ll and Klaus HornbergerS 
^Faculty of Physics, University of Duisburg - Essen, Lotharstrafle 1, 4^048 Duisburg 

We extend the theory of matter-wave interferometry of point-like particles to non-spherical objects 
by taking the orientational degrees of freedom into account. In particular, we derive the grating 
transformation operator, that maps the impinging state onto the outgoing state, for a general, 
orientation-dependent interaction potential between the grating and the molecule. The grating 
operator is then worked out for the diffraction of symmetric top molecules from standing light waves, 
and the resulting interference pattern is calculated in the near-held. This allows us to identify a 
signature of the orientational degrees of freedom in near-held matter-wave experiments. 
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I. INTRODUCTION 

Interference experiments with heavy molecules and 
nanoclusters amount to tests of the quantum superpo¬ 
sition principle at unprecedented scales w and may 
serve to measure molecular properties with high accu¬ 
racy m\- While it is an open question whether the 
superposition principle is valid at all scales [HHini , it is 
without doubt that the ro-vibrational degrees of freedom 
of increasingly large, non-spherical objects will at some 
point influence the interference signal. 

Far-held matter-wave experiments with heavy particles 
are challenging due to the small de Broglie wavelength, 
but near-held interferometry proved to be a powerful tool 
in the quest for high mass interference [HIIIIIHIS]. Near- 
held techniques are based on the Talbot effect, the repro¬ 
duction of the intensity pattern in the grating at certain 
distances further downstream. Since near-held interfer¬ 
ence effects are highly sensitive to even tiny phase mod- 
ihcations, it is natural to expect that the inhuence of 
the rotational state is most pronounced in the near held. 
In fact, signatures of the vibrational molecular dynamics 
have been observed in near-held experiments in 

In order to extend the established theory of matter- 
wave interferometry of spherical particles to non- 
spherical objects with orientational degrees of freedom, 
we draw on the results obtained for the dehection of ro¬ 
tating molecules [HHis]. There, it is a central result 
that the dehection angle of a rapidly rotating molecule 
is determined by its rotational state when entering the 
dehection held. In most cases of interest, the rotation 
of the molecule is initially thermally distributed, in par¬ 
ticular in the absence of a pre-aligning pulse. The dis¬ 
tribution of the rotational state then translates directly 
into a range of dehection angles, each observed with the 
thermal probability of the respective state. The distribu¬ 
tion describing the probability of a particular dehection 
angle is thus a central element in the theory of molecu¬ 
lar dehection and it will turn out to be similarly relevant 
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for the theory of matter-wave interferometry of rotating 
molecules. 

This article presents a full quantum theory of the in¬ 
huence of rotations in matter-wave interferometry, and 
we illustrate this theory by calculating the near-held 
interference pattern of symmetric top molecules in the 
Kapitza-Dirac-Talbot-Lau interferometer (KDTLI). We 
emphasize that the general formalism presented here is 
not restricted to the orientation state, but can be applied 
to other internal degrees of freedom as well. 

The article is structured as follows: In Sec. Ullwe de¬ 
rive the quantum mechanical grating transformations for 
the cases of a rotationally free and a rotationally dia- 
batic transit through the grating as well as their classical 
analogues. In Sec. |nT]we specify the grating transforma¬ 
tion for the dihraction of symmetric top molecules from 
a standing light wave in order to study the near-held 
interference pattern in the Kapitza-Dirac-Talbot-Lau in¬ 
terferometer (KDTLI), and we identify a signature of the 
rotational dynamics. We conclude in Sec. |IV| The ap¬ 
pendix provides a derivation of the distribution of classi¬ 
cal dehection angles of rotating symmetric molecules, as 
required in Sec. |III| 


II. THE GRATING TRANSFORMATION 

It is the aim of this section to determine the effect of 
the grating on the translational and rotational dynamics 
of the molecule. For this sake we consider a molecule of 
mass M that traverses with constant velocity Vz a dihrac¬ 
tion grating located at z = 0 and with the grating axis 
along the ^-direction. The assumption that Vz remains 
constant throughout the diffraction process is well justi- 
hed because in the experimental realization the kinetic 
energy of the longitudinal motion exceeds the average 
interaction potential and the kinetic energy of the trans¬ 
verse motion by orders of magnitude [50]. In addition, 
the restriction to a single velocity Vz is not a limitation 
because a hnite longitudinal coherence can always be in¬ 
corporated by averaging over the distribution of Vz in the 
end 

The orientation-dependent interaction between the 
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molecule and the grating is described by the grating po¬ 
tential V{x^z,Vl)^ where denotes the set of orienta¬ 
tional degrees of freedom (DOFs) of the molecule, such as 
the Euler angles, z = is the center-of-mass (CM) po¬ 
sition of the molecule in the flight direction at time t and 
X is the CM position in the grating direction. Since the 
extension of the particle beam in y-direction is usually 
small compared to the extension of the grating, we can 
neglect the y-dependence of the potential V{x,Vzt,^). 

In what follows, we distinguish between the rotation- 
ally free transit, where the (expectation value of the) 
rotational period r^ot of the molecule is constant dur¬ 
ing the passage and much smaller than the transit time 
TcM through the grating, and the rotationally diabatic 
scenario, in which the characteristic time tcm is much 
smaller than the rotational period T^ot- While the rota¬ 
tionally free transit is realized, for instance, in near-field 
matter-wave interferometry with laser gratings |12j . the 
diabatic transit can occur in far-held experiments with 
very thin material gratings m- 


A. Quantum Mechanical Description 

The central tool in the theory of matter-wave interfer¬ 
ometry of spherically symmetric particles is the grating 
transformation operator t, that maps the incoming trans¬ 
verse state p onto the outgoing transverse state p' = tpt^ 
with {x I p' \x') =t{x)t*{x'){x\p\x') [lailSlEO]. In the 
case of extended, non-spherical molecules, the grating 
transformation operator t must be adapted in order to 
account for the effect of the orientational DOFs. Let us 
start by deriving this operator from the time-dependent 
Schrodinger equation. 

The total Hamiltonian H contains, in addition to the 
grating potential V{x, z, Q), the CM kinetic energy Hcm 
in transverse direction x, as well as the free rotational 
Hamiltonian iLrot, whose form is determined by the sym¬ 
metry of the molecule |Hj. In what follows we keep the 
discussion general and denote by \£'m) the eigenfunctions 
of the rotational Hamiltonian with energy e^, where i la¬ 
bels the energy levels and m labels the degenerate states 
for each £. The symmetry of the rotor is arbitrary and 
so is the degeneracy vi for each 1. 

For the sake of a more compact notation we introduce 
for each i the tuple ^^(H) containing all the states with 
energy i.e. = (Dlfm). The length of the tuple 

(^(Q) is thus equal to ve. With these tuples the total 
wave function v['(x, Q,t) = (x, f2|'I'(t)) can be expanded 
as 


^'(x,Q,t) = ^e “^'^‘/^^(x.t)-^(Q). (1) 

e 


where x^(x,f) are the tuples of expansion coefficients 
iXg)m = Xim- Again, the length of X^(x,f) is 
Inserting the expansion Q into the time dependent 
Schrodinger equation with the total Hamiltonian H gives 


the coupled equations 


'ifidtXi,{x,t) = HcMXp{x,t) 


E' 


—lA 




Here, = et> — ee is the rotational energy-level 

spacing and we defined the grating potential matrix 


{Y,u,)mni'{x,Vzt) = Um V{x,Vzt,fl) i'm'j 


mension i/p x vgi. In addition, we denote the initial condi¬ 
tions to the Schrodinger equation Q by variables with¬ 
out the time argument, such as XpW- 

Equation ([^ describes the coupled time evolution of 
the expansion coefficients Ximix,t) due to the effectively 
time-dependent interaction between the molecule and the 
grating. The exact grating transformation t is given by 
the unitary time evolution of the system ([^; however, 
in practice a semiclassical approach is sufficient due to 
the small de Broglie wavelength in matter-wave exper¬ 
iments with heavy molecules |20j . In most cases it is 
even sufficient to determine the grating transformation t 
in the eikonal approximation, which can be regarded as 
the high energy limit of the semiclassical propagator . 
Physically speaking, the eikonal approximation treats the 
particle trajectories appearing in the semiclassical prop¬ 
agator as straight lines [22] • We now specify t explicitly 
for the rotationally free and for the rotationally diabatic 
transit through the grating. 


Free Rotor Transit 


In a rotationally free transit through the grating the 
molecule rotates rapidly during the passage through the 
grating, Trot ^ tcm, and the rotational energy clearly 
exceeds the average potential. Then the transverse CM 
wavefunctions Xi are nearly constant during the tran¬ 
sit and one can neglect the non-resonant terms in the 
Schrodinger equation ([^ (rotating wave approximation), 


indtXp{x,t)= HcM + Vp,{x,Vzt) X,ix,t). (3) 


The interaction potential is effectively diagonal in the 
angular momentum quantum numbers £ due to the fast 
molecular rotations. The corresponding expansion co¬ 
efficients ^(x,t) for different energies ee are mutually 
independent; however, in general the entries within each 
Xp{x,t) are coupled via Eq. @. 

In the eikonal approximation [111201 HD the scat¬ 
tered state ;^(x) behind the grating is according to the 
Schrodinger equation (1^ related to the impinging state 
Xpix) by x'(x) = '^{xy^ix) where the grating transfor¬ 
mation matrix is given by 


Z^{x) = exp 


i 


(4) 


It is a square matrix of dimension vp. The corresponding 
grating transformation operator t can be expressed in 
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terms of the matrix elements t™"* = {T^)mm' of the 
grating transformation matrix Q as 

(x) (g) (5) 

£ mm' 


Here, we included the aperture function \t{x, H)| G {0,1} 
describing the grating structure. Pure phase-gratings are 
characterized by |t(a;,H)| = 1 while for ideal interaction- 
free gratings the eikonal phase vanishes m- 

The grating transformation ([^ is valid for any interac¬ 
tion potential as long as the rotating wave approximation 
is justified. This is the case in most matter-wave experi¬ 
ments with large particles, such as near-held interference 
experiments with laser gratings dlHIS] or far-held exper¬ 
iments with thick gratings In Sec. Ill we will specify 
the grating transformation (O for laser gratings. In ad¬ 
dition, we note that in many practical cases the initial 
rotational state is a thermal mixture of angular momen¬ 
tum states and, hence, the grating transformation (§ can 
be regarded as the thermal average of angular momentum 
dependent grating transformations t™"* (x). 


Diabatic Transit 

In a rotationally diabatic grating passage the interac¬ 
tion time TcM is much smaller than the rotational pe¬ 
riod Trot and, in a classical picture, the orientation 
remains constant during the transit. In this case the 
quantum mechanical grating transformation gets diago¬ 
nal in the orientational coordinates H. In order to see this 
we note that for short transit times the rotating phases 
Au't/h ~ 2'KtlTT:ot in the Schrodinger equation ([^ can 
be neglected. This yields the coupled equations 


As anticipated, the diabatic grating transformation 
is diagonal in the orientational DOFs D. This coincides 
with the classical perception that the orientation is con¬ 
served during the diabatic passage through the grating, 
and is in contrast to the free rotor case ([^, which is di¬ 
agonal in the angular momentum quantum numbers 1. 
The diabatic transformation 0 can be appropriate to 
describe the transit through ultra-thin material gratings, 
e.g. made of graphene. 


B. Classical Grating Transformation 


In order to identify genuine quantum effects in near- 
held matter-wave interferometry, it is necessary to com¬ 
pare the quantum interference signal to the classically 
expected shadow pattern of the grating [T^] . This moire- 
type signal can in principle be obtained by solving Hamil¬ 
ton’s equations of motion for a rotating molecule in the 
grating potential V{x,v^t,Q). However, the problem is 
signihcantly simplihed by the classical analogue of the 
eikonal approximation. 

The classical state of the rotating molecule approach¬ 
ing the grating is described by its phase space distribu¬ 
tion function f{x,Px, H,pn), where pn denotes the vector 
of conjugate momenta to the angles fl. We seek the classi¬ 
cal grating transformation, that maps the incoming state 
f{x,px,^,pn) onto the outgoing state f'{x,px,^,pn) 
[ElES], i.e. the classical analogue of the grating trans¬ 
formation operator t mapping p onto p'. We discuss the 
free rotor scenario hrst. 


Free Rotor Transit 


ihdtX^{x,t) = HcMXf,{x,t) + E Kit 

I' 

( 6 ) 

which depend on the orientational DOFs only in para¬ 
metric fashion, since these equations are independent of 
the rotational energy levels £^. Defining the wavefunc- 
tion, 

='^)^{x,t) ■ (7) 


For a rapidly rotating molecule, the CM motion is de¬ 
termined by the grating potential averaged over a rota¬ 
tional period |24j . The resulting eikonal momentum kick 
Apx experienced by the molecule while passing through 
the grating reads as 

Apx{x,V,,pq) = -/ dt' / dz 

XrotVz Jo J — oo 

xdxV [x,z,n{t')]. (10) 


allows us to rewrite the coupled Eqs. Q in the form 


ihdt^{x, O, t) 


Hcm + V (x,Vzt,fi) 


<^(x,fi,t), 


( 8 ) 


which can now be solved in the eikonal approximation. 
The diabatic grating transformation operator t, mapping 
the initial state H) onto the outgoing state H), 
can now be written in the eikonal approximation as 


i = \t{x, H)| exp 


i 

Hvz 


dz V{x, z, Cl) 


(9) 


The transferred momentum (101 is a function of both the 
transverse CM coordinate x and of the initial orientation 
state (H,pf 2 )) that determines the rotational dynamics. 
In addition, it is reasonable to neglect the influence of 
the grating on the rotational dynamics, since the rota¬ 
tional energy is much higher than the average interac¬ 
tion potential. The free rotor approximation is in most 
practical cases well justified due to the high rotational 
temperature of the molecules in the experiments [5]. 


With the help of the eikonal momentum kick ( |Io| , we 
can express the outgoing distribution f'{x,Px, ^,Pn) as a 
CM momentum convolution of the impinging distribution 
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f{x,p^,n,pn), 


f{x,p^,n,pn) = j dp'^ Tci{x,p^ -p'^,n,pQ), 


xf{x,p'^,n,pn), 


( 11 ) 


with the grating transformation function 

Tci {x,poo,fl,pn) = |t(a:, [px - Apx{x,fl,pn)] (12) 

The function Tc\{x,Px,^,Pn) is the classical analogue of 
the quantum mechanical grating transformation opera¬ 
tor t. We remark that the transformation (11) conserves 


the angular momenta pn, as it was the case in the quan¬ 
tum mechanical case We will now specify the grat¬ 
ing transformation function for the rotationally diabatic 


case. 


Diabatic Transit 


In the case of a rotationally diabatic transit, the ro¬ 
tational period is much longer than the transit time and 
thus the orientation of the molecule can be considered 
as being constant during the passage through the grating. 
The resulting CM momentum kick (101 is 


1 f°° 

Apx{x,il.) = -/ dz dxV {x,z,il.), 

'^Z J — oo 


(13) 


which is a function of the transverse CM position x and 
of the orientation Q. In a similar fashion, we obtain the 
eikonal angular momentum kick ApQ, 

Apn = -—[ dz dnV {x,z,n), (14) 

Vz 7-00 


which is also a function of x and 17. The resulting 
grating transformation is now given by a CM momen¬ 
tum and an angular momentum convolution analogous 
to Eq. © with the grating transformation function 
Tci{x,px,^,pn), 

Tc\ {x,Px,^,Pq) = |t(x, 17 ) 1(5 [px - Ap2,(a::,f7)] 

X(5[pn - Apo(a;, 17)]. (15) 

Having derived the quantum and classical grating trans¬ 
formations, we can next apply them to the molecular 
diffraction from standing wave laser gratings. 


III. INTERFERENCE OF SYMMETRIC TOP 
MOLECULES IN THE KDTLI 

Here, we first discuss the diffraction of rapidly rotating 
symmetric top molecules from laser gratings in order to 
illustrate the previously derived grating transformation. 
The obtained transformation operator is then used to de¬ 
termine the quantum fringe visibility as well as the classi¬ 
cal shadow pattern in the Kapitza-Dirac-Talbot-Lau in¬ 
terferometer (KDTLI) [3 [TTl [T2] . The KDTLI is a near- 
held interferometer consisting of three gratings that all 


share the same grating period d. The Hrst and third grat¬ 
ing are material masks, while the central one is a stand¬ 
ing light wave. The transverse coherence of the incoming 
particle beam is prepared by the hrst material grating 
at distance L in front of the standing wave. Diffraction 
occurs at the central grating and the signal is detected 
with the help of the third one at distance L further down¬ 
stream. The KDTLI is currently the working-horse for 
high mass interference experiments in Vienna [2]. 

For the theoretical description of rotating molecules in 
the KDTLI it is a reasonable approximation to neglect 
the inhuence of the hrst and the third grating on the 
orientational DOFs because the diffraction relevant for 
interference takes place only at the central grating. The 
total transverse state operator p behind the third grat¬ 
ing can be obtained by successively applying the grating 
transformations of the three gratings as well as the in¬ 
termediate unitary evolutions [niEo]. Finally, the ori¬ 
entational DOFs are traced out in order to obtain the 
interference pattern on the screen. 


A. Standing Wave Grating Transformation 


We now evaluate the grating transformation t for sym¬ 
metric top molecules traversing a standing-wave laser 
grating. This operator will then be used to calculate the 
near-held interference pattern of symmetric molecules in 
the KDTLI, but it can also be applied to other situations. 
We consider a polarizable, symmetric molecule with mo¬ 
ments of inertia I = Ii = I 2 and I 3 {I /13 > 1 for prolate 
particles and 1/2 < I/I 3 < 1 for oblate discs), that is 
diffracted from a Gaussian standing laser wave of wave¬ 
length A. The laser wave is linearly polarized in the di¬ 
rection n and acts as a pure phase grating, \t{x, I7)| = 1. 
The intensity of the Gaussian standing laser beam aver¬ 
aged over one optical cycle is given by [12] 


Iix,z) = 


8P 

TnUyWz 


/ 2z2^ 

\ ■ 2 ( ^\ 


sm ( TT— 

V wi j 

1 V dJ 


( 16 ) 


where d = A/2 is the grating period, P the laser power 
and Wz is the beam waist in z-direction. Since the ex¬ 
tension of the incoming particle beam in y-direction is 
usually small compared to the beam waist in y-direction 
[TT] it is natural to neglect the w-dependence of the in¬ 
tensity (@. 

Denoting by ay and a_L the two independent compo¬ 
nents of the polarizability tensor of the particle (along 
its symmetry axis and perpendicular to it, respectively), 
the grating potential can be expressed as dlillHISI Eij 


V (x,z,d) = - 


4P 


■ exp 




neoCWzWy 
X (ay — Aofsin^d) sin^ ' 


( 17 ) 


Here Aa = a\\ — a_L is the polarizability anisotropy of the 
molecule (Aof > 0 for prolate particles) and 9 is the nu¬ 
tation angle with respect to the field polarization n. The 
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grating transformation can be safely determined in the 
free rotor approximation, Eq. ([^, since the laser beam 
waist along the flight direction is approximately Wz = 20 
/rm mmi]. For an exemplary molecule (diazobenzene 
with mass M ~ 1030 amu and length Lmoi — 3.5 nm) 
with velocity Vz = 100 m s“^ and rotational temperature 
T = 600 K, we thus have Trot/TcM ~ 10“^. 

In order to calculate the grating transformation ma¬ 
trix Q we need to specify the grating potential ma¬ 
trices, whose particular form depends on the symme¬ 
try of the rotor. Here we restrict our discussion to 
symmetric top molecules for reasons of simplicity. The 
eigenstates \imk) of the symmetric rotor (with classical 
Hamilton function (Al), see appendix) are labeled by 
the three quantum numbers £ € Nq, m € 
and k S {—£, ■■■,£} with eigenenergies e^k [Ml 127] ■ De¬ 
noting by € [0,27r), 6 S [0,7r) and ip G [0,27r) the 
three Euler angles with respect to the field polarization 
n {z-y'-z" convention), the configuration space represen¬ 
tation of the states \imk) can be given explicitly in terms 
of the Wigner H-matrices. In particular, the eigenfunc¬ 
tions are {^\^mk) = \/2ATT/(27r-\/2)Il^^(H) [Ml HZ], 
where is related to the (small) Wigner d-matrix 

by We remark that the 

(small) Wigner d-matrices d^j. (d) are real due to the em¬ 
ployed definition of the Euler angles [22] • The basis kets 
\£mk) are complete and orthonormal with respect to the 
infinitesimal volume element dfl = dip dd dip sin d. 

The interaction potential 0 is a function of the nu¬ 
tation angle d only and thus the quantum numbers m 
and k are conserved. The resulting grating potential is 
diagonal in all three quantum numbers £, m, fc, and the 
diagonal elements can be given with the help of the expec- 


sin^ d 


imk ). Expressing 


tation values Qtmk '■= {imk 
these expectation values in terms of Wigner id-matrices 
[m and using the properties of the Wigner 3-j symbol 
[23127] yields the algebraic form 

^ 1,1 (2mf + {2kf - 1 

Qirak - 2 + 2 ~l^r- 1) (2T+‘Sj" 

3 {2mkf 

2i{i+l){2i-l){2e + 3)' ^ ’ 

We remark that in the limit of a linear rotor, ///a —^ oo 
and thus k = 0, the well-known expectation value |14] 

^ _1 1 ( 2 to ) 2-1 

- 2 + 2 ( 2 £_ 1 )( 2 ^ + 3 )’ 

of the linear rigid rotor is recovered. 

We are now ready to identify the grating transforma¬ 
tion operator for symmetric top molecules traversing a 
standing wave laser grating by inserting the expectation 
value (181 into the grating transformation matrix ([^. 
The resulting operator ([^ is 


e 


t = 




^ 0 exp 

'^\imk) {£mk\, 


1 - 


Aa 


•mk sm 


X 

’"d 


where (po = 4:a\\P/eochwyVz\/^ is the eikonal phase [T^] 
defined with the polarizability ay along the symmetry 
axis of the molecule. It is important to note that the 
eikonal phase imprinted on a particle during the grating 
passage depends on all of its angular momentum quan¬ 
tum numbers £, m and k. The final signal, obtained by a 
trace over the orientational DOFs, is thus an average over 
signals from different grating transformations (20), each 


weighted with the probability of the corresponding angu¬ 
lar momentum state \£mk). This matches with the fact 
that the classical deflection angle of molecules traversing 
an electrostatic field depends on the angular momentum 
of the deflected particle [HHis]. For nearly isotropic par¬ 
ticles, Aa/ay 1, the transformations (20) are all equal 
and the average over angular momentum states can be 
neglected. 


B. Quantum and Classical Fringe Visibility 


Having specified the grating transformation (20), the 


quantum fringe visibility V of symmetric top molecules 
in the KDTLI can be calculated by applying the trans¬ 
formation ([^ with ( [20| to the quantum phase space for¬ 
malism presented in |12j . The common period of all three 
gratings is denoted by d and the de Broglie wavelength of 
the incoming rod by Ads = h/Mvz- The Talbot length, 
the characteristic length scale in near-held interferome¬ 
try m, is Lt = d^/AdB- A straightforward calculation 
yields the sinusoidal quantum fringe visibility 


V = 2sinc^(7r/)y^ ^ ptmk 

i—0 m,k——£ 


X J 2 


1 - 


Aa 

- ( 

ay 


fc sm TT 


Tx 


, ( 21 ) 


( 20 ) 


where / is the opening fraction of the Hrst and the third 
grating, J 2 (-) is the second order Bessel function of the 
hrst kind, L is the distance between the gratings and 
Pimk is the statistical weight of the angular momentum 
state \£mk). The interference contrast can be regarded 
as the average of point particle visibilities with (£, m, k)- 
dependent eikonal phases and weights pimk ■ 

Since the molecules are emitted from a thermal source 
into vacuum the rotational DOFs follow a thermal distri¬ 
bution, pimk ~ exp(—e^fe/ZcBT), at a very high tempera¬ 
ture, k^T ^ h?/1. Then the sum over angular momenta 
in Eq. ( [^ can be replaced by the integral over the 
corresponding classical distribution [27] and Eq. (21) is 
further simplihed. In particular, we denote by pth( 9 ythe 
probability density of the variable q = Q{Ej:ot,Pip,Pip), 
where Q{E^ot,Ptp,Pi!) is the classical free temporal mean 
value of sin^ 9{t) depending on the conserved rotational 
energy E^ot and on the canonical momenta Py, and Pif, 
of p> and Ip rotations, respectively. A simple expression 
for Q(£'rot,P(^,Pi/,), as well as for the thermal distribu¬ 
tion Pth ( 9 ), is derived in the appendix. This probability 
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FIG. 1. (Color online) The thermal distribntion Pth(9), Eq. 
\22\ , of the temporal average q of sin^ 9{t) for the free sym¬ 
metric rotor with different moments of inertia, I/I 3 = 1/2 
(solid line), I/I3 = 10 (dashed line) and I/I3 —>■ 00 (dot- 
dashed line). 


density is depicted in Fig. and reads 


Pthiq) = \li^l du 

■ils Jciq) 


1-(1--,.- 

- 1/2 


1 - 3/2 


--A{u^ + l- 2 q) 


( 22 ) 


where the integral must be taken over the union of two 
intervals, ({q) = Ci{q) U The first interval Ci{q) = 

[0, ^min[A(g)]], where A{q) = {1/3, 2q — 1,1 — q}- This 
contribution to the distribution (221 vanishes for q < 


1/2. The second interval ( 2 {q) = [-\/max[^(( 7 )], 1]. The 
distribution ( p^ depends on the fraction I/I3 only and 
is independent of the thermal energy k^T. 

In Fig. we show the distribution (22) of the sym¬ 
metric rotor for the oblate limit {I/I3 = 1/2), a prolate 
particle {I/I3 = 10) and the linear rotor {I/I3 —>■ 00 ). 
Similar figures were obtained numerically in |15j . For 
finite I/I3 the probability density (22) is discontinuous 
at g = 1/2, which follows from the definition of the set 
C(g) and it diverges at g = 2/3, as can be observed di¬ 
rectly from Eq. (22). In the limit of the linear rotor, 
I/I3 —)■ 00 , the established [H] form pth(g) = l/V^g — 1 
is recovered. 


Using the probability density (22) the quantum fringe 


visibility V takes on its final form 

V = 2sinc^(7r/) / dg pth(g) 
Jo 


X J2 


0 . I 1 - dn A) 


(23) 


niently calculated with the help of the phase space grat¬ 
ing transformations (11). The classical momentum kick 


(10) transferred to the molecule by the grating potential 
(17) is 


Ap,c(x,Erot,Pq,) = 


7rh(j)Q 


Arv 

1 Pipt P'lp) 

ail 


X sin 


(-1) ■ 


(24) 


Following the treatment in 
fringe visibility 


one obtains the classical 


Vci = 2sinc^(7r/) / dg pth(g) 
Jo 


X J2 


Aa \ 

1- q —r — 

a\\ J Lt 


(25) 


where we assumed the orientational DOFs to be ther¬ 
mally distributed. 

The classical visibility (25) decays as with in¬ 

creasing grating separation L (and decreasing particle 
velocity Vz) while the quantum visibility ( |23[ ) is periodic 
in L/Lt [mUl]. This can be used to discriminate be¬ 
tween genuine quantum behavior and classical shadow 
effects m as illustrated in Fig. for the linear rotor 
///a —>■ 00 . In Figs. and we consider an exemplary 
molecule (M = 1030 amu, Lmoi = 3.5 nm, a = 47r£o x 50 
Vz = 100 ms I/I3 ~ 00 [11] ) traversing the 
KDTLI {d = 266 nm, L/Lt = 0.5, / = 0.42, Wz = 20 pm 
[HHH]) for three different relative anisotropies Aa/ay. 

An experimentally observable signature of the orien¬ 
tational DOFs can be found in the absolute value of the 


quantum fringe visibility (23) as a function of laser power 
P. This is illustrated in Fig. ^ While the value of subse¬ 
quent maxima in the visibility are strictly decreasing for 
spherical particles, this is not the case for non-spherical 
molecules. For large relative anisotropies Aa/ay the 
thermal average over angular momentum states leads to 
the appearance of additional side peak between the ma¬ 
jor recurrences. This is a signature of the orientational 
DOFs of the molecule traversing the grating. 

In Fig. 1^ we show the absolute value of the quan¬ 


tum fringe visibility (23) as a function of laser power 


for three differently shaped prolate molecules. All other 
parameters are as for the linear molecules of Fig. 
(Aa/ay = 0.9). While the visibility coincides with the 
visibility of the linear molecule for large ratios I/I3, it 
approaches the behavior of spherically symmetric object 
for I/I3 —)■ 1. The signature of the orientational DOFs 
discussed above is most pronounced for linear molecules 
but can be observed for all prolately shaped molecules. 


IV. CONCLUSION 


In order to identify genuine quantum interference effects 
we must compare the visibility ( |23[ ) to the visibility Vd 
of the classical shadow pattern |12j. which is most conve- 


We extended the theory of matter-wave interferometry 
to large, non-spherical particles by accounting for the 
influence of the rotational dynamics. In particular, we 
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FIG. 2. (Color online) Quantnm (upper curves) and classical 
(lower curves) absolute sinusoidal fringe visibility as a func¬ 
tion of relative separation L/Lt for a linear molecule in the 
KDTLI for three different values of Aa/a\\. 



FIG. 4. (Color online) Absolute value of the sinusoidal quan¬ 
tum fringe visibility ( |23| l as a function of laser power for the 
prolate symmetric top in the KDTLI for different values of 
I /13 (Aq/q|| = 0.9). 



ing, and we showed how it enters the description of sym¬ 
metric top particles in the Kapitza-Dirac-Talbot-Lau in¬ 
terferometer as performed at the University of Vienna 
m- In these experiments the typical transit time ex¬ 
ceeds the rotational period by orders of magnitude, mak¬ 
ing the transit rotationally free. A signature of the rota¬ 
tional dynamics was pointed out in the predicted quan¬ 
tum fringe visibility as a function of laser power. We also 
derived a closed-form expression for the distribution of 
deflection angles in classical deflection experiments with 
symmetric top molecules, as required in this context. 


V. ACKNOWLEDGMENTS 


FIG. 3. (Color online) Absolute value of the sinusoidal quan¬ 
tum fringe visibility (211 as a function of laser power for a 
linear rigid molecule in the KDTLI for three different values 
of Aa/ay. 
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derived the grating transformation operator for the ro¬ 
tationally free and for the rotationally diabatic transit. 
This operator describes the modihcation of the transverse 
quantum state due to the orientation-dependent interac¬ 
tion with the grating. In addition, the classical shadow 
pattern was derived in order to provide the tools required 
for the identihcation of genuine quantum effects in near- 
held matter-wave interferometry. 

If the molecule rotates rapidly with high energy, the 
transit is rotationally free and the grating transformation 
depends only on the angular momentum of the impinging 
particle. On the other hand, if the transit time is much 
shorter than the average rotational period, the grating 
transformation depends on the orientation of the particle. 

We worked out the grating transformation for symmet¬ 
ric top molecules traversing a standing-wave laser grat- 


Appendix A: Rotation Statistics of Symmetric Top 
Molecules 


In this appendix we present a derivation of the thermal 
probability distribution /th(?') of classical realizations r 
of the temporal mean value of cos^ 0 {t) for symmetric top 
molecules, as required in Sec. Ell This distribution is of 
central interest in the theory of molecular deflection ex¬ 
periments because it translates directly into the distribu¬ 
tion of deflection angles [HIIMS] and it has been eval¬ 
uated numerically with the help of Monte-Carlo methods 
in m- The distribution ( [2^ of realizations q of sin^ 9{t) 
can be trivially obtained by substituting r = 1 — 9 in 
fth{r). 

We follow the notation of Sec. IIIII and denote the mo¬ 
ments of inertia of the symmetric top hy Ii = I 2 = I 
and I 3 . The three orientational DOF D = {ip,9,'4’) are 



















the Euler angles in the z-y'-z" convention, and the re¬ 
spective conjugate momenta are pQ = (p^^pg^p^). The 
classical Hamilton function H,.ot{^,Pn) of the free sym¬ 
metric top is given by [53] 

= + + (Ai) 

zi \ sm 0 J 2 I 3 

and there are three conserved quantities: the total energy 
E^ot = H^oti^^Pa), as well as the angular momenta p^ 
and p^. Hence, the equations of motion are integrable 
and the period of 0-rotations acquires the form 


2ttI 


'^rot — 


^ 2 E,,J+pl{l-I/h)' 


(A2) 


with the partition function Z = IkBT^/2TTI 

[28], the distribution fth{r) of r = R{E^ot,Pip,Piij) is given 
by 

/th(r) = J dT S[r - R{E^ot,Pip,P^)]pthi^,Pn), (A7) 

where the integral covers the whole phase space (dT = 
dp dO dtp dp^ dpg dp^). On the ot her h and, one can also 
write the distribution fthir), Eq. (A7|, in terms of the 
probability density (/th(wi,U 2 ), 


qth{ui,U2) = J dT S (ui 


Pip'^rot \ 


27rl ) 

x<5 (ui - ^^1^) Pth(0,pn), (A8) 


according to 


In the limit of very prolate molecules, ///a —> 00 , "0- 
rotations do not contribute, = 0, and the period 
Trot approa ches the period of the linear rotor, Trot —>■ 
7r^2J/£;rot- 


Separating Hamilton’s equations with (Al) yields the 
temporal average of cos^ 0(t) over one rotational period 


\ r'^Tot 

R{E^ot,Pv,P 4 >) = - / cos^ 9{t') 

^rot Jo 


pI+pI 


1 1 

2 " 2TErot/+p|(l-///3) 


Pippljj 


2Erot/+p|(l-///3) 


(A3) 


We observe that the quantum mechanical expectation 
value Rimk = 1 — Qemk, Eq. (181, closely resembles this 
expression. 


fth{r) = J dui 


duo S 


/-I 


r-R{ui,U 2 ) qth{ui,U 2 ). 

(A9) 

The integral (A 81 can be evaluated by defining the di¬ 
mensionless quantities = {p^ —p-ip cos 9) / sin 9\/Ik^T, 
ujg = pgjy/ThsT and w,/, = Pijjl\/Ik^T and transforming 
the vector (uj^,u!g,uj^) to spherical coordinates, {u!,^,r]) 
(where ui = Prot/I)- This gives 


qth{ui,U 2 ) = y- 

4 V ^3 




ul 


-1 - 3/2 


,(A10) 


which is independent of the temperature T, uniformly 
distributed in ui, a function of the ratio / /I 3 only, and 
normalized by construction. 


Finally, inserting the distribution (AlO) into expres¬ 


sion (A9| yields the thermal distribution /th(?’ 


fth{r) = 


3/3 Jcir) 


du 


’-'’-ft)”’ 

- 1/2 


-1 - 3/2 


A more compact expression for the temporal average 
(A3) is obtained in terms of the relative frequencies of p- 

X 

-f 2 r - 1 ) 

and 0-rotations 


L V ■5/ J 


Pif'^rot J P4>'^Tot 


(A4) 


which satisfy —1 < ui,U 2 < 1. In particular, introducing 
the angular momentum scale p^ot = SttJ/ Tj-pt we observe 
that IPipy] < Prot- Then, the average ( ]A3| ) can be written 
as 

R{ui,u 2 ) = ^-^{ul + ul) + ^ulul, (A5) 

which will be of advantage in what follows. 

If the orientational DOFs are distributed in phase 
space according to the thermal distribution 


Pth(f^,Pn) = ^ exp [-i7rot(f^,Pn)/fcBT], (A6) 


(All) 


Here, the set C(r) C [0,1] denotes the fraction of the unit 
interval [0,1] where the integrand is real. In particular, 
for arbitrary r £ [ 0 , 1 ], C(r) = Ci(’’) U 02 (?') consists of 
two intervals Ci, 2 {r)' the first is Ci(^) = [0, y/min[A(r)]], 
where A(r) = {1/3,1 — 2r, r}. This contribution van¬ 
ishes for r > 1/2. The second interval C 2 (?') is C 2 ('c) = 
[-ymax[A(r)], 1] to one. The probability density function 
(22) of the temporal average q = (sin^ ^(0)i ^-s required 
in Sec. |HI[ is obtained by the substitution q = 1 — r. 


Finally we note that the distribution fthir) (All) ap¬ 
proaches in the rigid rotor limit / /I 3 —^ 00 (in agreement 

with [33]) 


/th(r) = 


1 




(A12) 
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